Solvent-induced micelle formation in a hydrophobic interaction model 
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We investigate the aggregation of amphiphilic molecules by adapting the two-state Muller-Lee- 
Graziano model for water, in which a solvent-induced hydrophobic interaction is included implicitly. 
We study the formation of various types of micelle as a function of the distribution of hydropho- 
bic regions at the molecular surface. Successive substitution of non-polar surfaces by polar ones 
demonstrates the influence of hydrophobicity on the upper and lower critical solution temperatures. 
Aggregates of lipid molecules, described by a refinement of the model in which a hydrophobic tail 
of variable length interacts with different numbers of water molecules, are stabilized as the length 
of the tail increases. We demonstrate that the essential features of micelle formation are primarily 
solvent-induced, and are explained within a model which focuses only on the alteration of water 
structure in the vicinity of the hydrophobic surface regions of amphiphiles in solution. 

PACS numbers: 



I. INTRODUCTION 

In an aqueous solution, the hydrophobic parts of am- 
phiphilic molecules tend to separate themselves from wa- 
ter molecules by forming aggregates, such as micelles 
and microcmulsion droplets. The simplest amphiphilic 
structure occurs if the polar and hydrophobic parts of 
the amphiphilic molecule are well separated into head 
and tail regions. Motivated by the structure of sphin- 
golipids and glycolipids (Fig. QJ, we will refer to this 
type of molecule as 'lipids' (although we note here that 
this generic structure is not common to all classes of 
lipid). For molecules such as those in Fig. ^ micelles 
consist of a polar outer surface and a hydrophobic core 
which contains all the tails. By contrast, amphiphilic 
species whose polar and hydrophobic regions are dis- 
tributed over the entire molecule, rather than being 
clearly separated, aggregate to form more general assem- 
blies for the purpose of minimizing the hydrophobic area 
per molecule exposed to the aqueous phase. We will re- 
fer to this general category of mixed HP molecules as 
'amphiphiles.' Poly(N-isopropylacrylamide), a polymer 
belonging to this class, exhibits a phase transition at a 
lower critical solution temperature (LCST) from a ho- 
mogeneous solution, where the polymers are completely 
soluble, to a system of two separated phases pj. At an 
upper critical solution temperature (UCST) the organic 
phase disaggregates, and above this temperature the am- 
phiphilic molecules are again soluble due to entropic ef- 
fects |2| . Substitution of polar by hydrophobic monomers 
in amphiphiles of given length leads to alterations of the 
critical solution temperatures which depend on the size 
of hydrophobic surface regions . 

Micelle formation, and the aggregation of amphiphilic 
molecules in general, may be treated as a phase separa- 
tion occurring at a critical micellar concentration (CMC) 
which describes the density of amphiphiles where the sys- 
tem enters the two-phase region Above the CMC, 
amphiphilic molecules in aqueous solutions self-associate, 



forming small aggregates to decrease the net contact be- 
tween their hydrophobic surfaces and the surrounding 
solvent. A small fraction of the lipid molecules remains 
free in the solution, with a concentration close to the 
CMC value. As the concentration increases, the onset of 
a semi-dilute regime is found, where the system may be 
considered as a solution of relatively few water molecules 
dissolved in an amphiphilic medium p|. 

For lipid molecules, the CMC decreases as the length 
of the hydrophobic chain increases [|| 0, H, @ , indicat- 
ing a stabilization of aggregates as a consequence of the 
stronger net repulsion between the tail and the surround- 
ing water. For the same reason, the LCST is thought 
to decrease as the tail grows, as suggested in Ref. |2j. 
One means of probing the nature of the phase diagram 




FIG. 1: Space-filling models of typical lipid molecules: (a) a 
sphingolipid (sphingosine, C1SH37O2N), with one non-polar 
hydrocarbon tail; (b) a glycolipid (C40-HV2O8-/V), with two 
tails. The relatively large, polar head-groups contain in addi- 
tion oxygen and nitrogen atoms. 
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and the effective hydrophobic interactions would be by 
systematic alteration of the polarity of the amphiphilic 
polymers in solution. Many theoretical and experimental 
studies have been conducted to describe the various ag- 
gregation phases of amphiphilic molecules in aqueous so- 
lutions 

However, no comprehensive investigations have yet been 
performed concerning the influence of the distribution 
of polar groups in amphiphiles on the aggregation phase 
diagram, or concerning the mechanism underlying the 
process of self-aggregation. 

The aim of this study is to substantiate experimental 
results indicating a decrease in LCST as the degree of 
hydrophobicity increases, and to analyze the dependence 
on density and hydrophobicity of the aggregation phase 
diagram. We will conclude that the principal properties 
of amphiphiles in aqueous solutions are solvent-induced, 
in that they are explained by alterations in the struc- 
ture of liquid water in the vicinity of the hydrophobic 
regions of solute particles. We begin by introducing a 
simple hydrophobic-polar (HP) description (Sec. Ill A|l of 
amphiphilic solute particles (Sec. Ill B|) on a cubic lat- 
tice, and then extend the model to describe lipid par- 
ticles of varying tail length (Sec. Ill C|l . In Sec. IIII Al 
we investigate the changes in the phase diagram asso- 
ciated with an increasing proportion of polar groups in 
amphiphilic molecules and with changes in their distri- 
bution. Sec. IIII 61 presents a similar analysis for the for- 
mation of lipid micelles in aqueous solutions. In Sec. IIVI 
we discuss the implications of our results and provide a 
brief conclusion. 



II. MODEL AND METHODS 

A. Water Structure 

The solvation of amphiphiles in aqueous solutions and 
their self-association into micellar aggregates are gen- 
really considered as a consequence of the effective hy- 
drophobic interaction between polar water and the non- 
polar regions of the solute molecules |19| . The unique 
properties of the aqueous medium which generate this in- 
teraction arise from the ability of water molecules to form 
strong hydrogen bonds both among themselves and with 
the polar groups of solute molecules. The formation and 
disruption of extended hydrogen-bonded networks leads 
directly to the delicate balance of enthalpic and entropic 
contributions to the solvent free energy which is respon- 
sible for the existence of a closed-loop aggregation regime 
in many solutions of polar molecules. This behavior is en- 
capsulated in an adapted version of the model of Muller, 
Lee and Graziano (MLG) |20l |21| . whose key features 
are as follows: consideration of a coarse-grained system 
whose sites contain clusters of water molecules of a size 
which matches that of the solute particles; a bimodal dis- 
tribution of water clusters reflecting "ordered" sites with 
mostly intact hydrogen bonds, and "disordered" sites 




FIG. 2: Energy levels of a water site in the bimodal MLG 
model. The states are denoted os — ordered shell (cage con- 
formation), ds — disordered shell, ob — ordered bulk, and db 
= disordered bulk. 



with relatively fewer intact hydrogen bonds; a further 
subdivision into "shell" sites neighboring the solute par- 
ticles and "bulk" sites which are not disrupted by prox- 
imity to a polar surface; a set of microscopic energy and 
degeneracy parameters for these sites whose sequences 
are determined by the fact that breaking of hydrogen 
bonds increases site enthalpy but simultaneously raises 
the site degeneracy (number of orientational degrees of 
freedom of the water molecules), as represented in Fig. EI 
Full details of the foundations and qualitative properties 
of the model can be found in Refs. [23] and [23|. 

The microscopic origin of the energy and degeneracy 
parameters for the four different types of water site (or- 
dered/disordered and shell/bulk) in the adapted MLG 
model, based on both experimental |2J, |2|| |26j and the- 
oretic al [27Lp3. l29l l30| analysis, is discussed in detail in 
Refs. [22|,[31j. Here we highlight the competition of en- 
thalpy and entropy terms using only the example of the 
ordered states. Although the insertion of a hydropho- 
bic molecule leads to a destruction of local hydrogen- 
bonding, at low temperatures the water molecules are 
found to rearrange in a cage-like structure around the 
solute molecule, which because of the orientational ef- 
fect of the polar surface is formed by stronger hydrogen 
bonds. The consequent net ener gy reduction results in 
dissolution of the solute particle |24l I2H l2r| . At higher 
temperatures, however, the additional entropic contribu- 
tions available from the bulk solvent favor a minimization 
of local water restructuring, which drives the aggregation 
of hydrophobic solute particles to minimize their total 
surface exposed to water. 

While it is possible to find perfect cages of hydrogen 
bonds for small hydrophobic particles, the formation of 
complete cages around large solute particles is prevented 
sterically. However, such large solute particles, and also 
their aggregates, may instead be surrounded by a number 
of partial cages, depending on their shape and surface 
roughness |2g. Atomic groups exposed at the surfaces 
of large molecules and aggregates present locally curved 
surface structures |32( which allow for a confined forma- 
tion of partial cages. The lipid molecules we analyze (see 
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Fig. 1) have a characteristic lengthscale of order 0.3 nm 
for the diameter of the chain and any side groups, 0.5 
nm for the polar head, and a maximum length of 2-3 
nm for the longer chains. The chain and head diame- 
ters are the lengthscales which determine the curvature 
of the molecule, and thus the possibility for cage for- 
mation. The considerations of the previous paragraph, 
specifically low-temperature solution, are then expected 
to remain valid for amphiphilic and lipid molecules, and 
this is confirmed by the success of models based on these 
concepts in describing large molecular species |33ll33.l35| . 
That the analysis is applicable for solute particles with 
a certain degree of curvature and surface roughness does 
not appear to be a significant restriction for molecules in 
the size range of most interest for amphiphilic and lipid 
characteristics. 

Returning to the parameters of Fig. [21 we stress that 
the properties of the sysytem are critically dependent 
on the sequences E ds > E d b > E Q b > E os of energy 
levels and q ds > qdb > q Q b > q os of the correspond- 
ing degeneracies, which describe the entropy. These se- 
quences are confirmed by a range of experimental mea- 
surements, and indeed if they are not maintained the so- 
lution does not exhibit a closed-loop aggregation regime. 
However, the results of the calculations to follow are not 
particularly sensitive to the exact values of the differ- 
ences between these parameters, and are of course in- 
dependent of their absolute values. We have used the 
energy values E ds = 1.8, E db = 1.0, E ob = -1.0, 
and E os — —2.0, which are thought to be qualita- 
tively representative for aqueous solutions, and which 
have been successful in describing different types of so- 
lution H3, HH, |3f| . The respective degeneracies, nor- 
malized to a non-degenerate ordered shell conformation, 
are taken to be q ds = 49, q c ib — 40, q b — 10 and 
q os = 1 H3- These relative values have been found 
to be appropriate for reproducing the phenomenology 
of hydrophobic interactions p2l , protein denatura- 
tion [23, |34( , swelling of biopolymers j3f| , and cosolvent 
effects on solubility of hydrophobic particles [23, 36] . Pre- 
cise values of the microscopic parameters may in fact be 
refined by comparison with experimental measurements 
to yield semi-quantitative agreement for different solu- 
tions jU 03, E^, |3(| . The energy scale is correlated di- 
rectly with a relative temperature scale, which we define 
as k B T ee (3- 1 . 



B. Amphiphiles 

To include not only purely hydrophobic solute par- 
ticles but also amphiphilic molecules with varying 
conformations of polar and non-polar regions, as in 
the experiments of Ref. 0, we represent the particles 
as cubes on which each face may be either polar or 
hydrophobic. A neighboring water site, which is homoge- 
neously polar, interacts only with the side of the particle 
which is oriented in its direction. If this side represents 



a polar group, the neighboring water site is in a bulk 
state, whereas if the side is hydrophobic the water site 
is considered a shell state. Thus polar sides and water 
are considered as having the same effect on neighboring 
water molecules. However, polar faces do not contribute 
to the free energy. In this coarse-grained model, one side 
of a site may represent more than one chemical group, 
and thus corresponds to a net characterization of the 
surface area of the solute molecule under consideration. 

Hydrophobic solute particles are generally larger than 
water, and thus a water site in the model consists of a 
group of molecules. On a lattice where each site has z 
nearest neighbors, the energy of a system of N sites, occu- 
pied either by amphiphiles (n^ = 0) or by water (rii = 1), 
is given by the Potts-like Hamiltonian 



H[{n i },{a i } ) {k}} = 



N 



(E os Si :<To3 +-Erf s 5j |tTd jy^(l- Xi^(k) ) 



(E ob 5i 



Edb8i,<j db ) ^2 A 

(ji) 



(1) 



where \ : j(k) is a side variable depending on each of the 
nearest neighbors j of site i, and takes the value 1 if the 
neighboring side of site j is a polar side of a solute parti- 
cle or water, and otherwise. The variable k takes values 
from 1 to n c j, where n c j enumerates all equivalent ori- 
entations of site j. On a cubic lattice, the total number of 
sides of a solute particle is z — 6 = Np + Njj, where Njj 
is the number of hydrophobic cube faces. Thus for the 
cubic lattice n c j = 1 for water and for solute particles 
with Njj — 6, n c j — 6 for Nh — 5, n c j — 3 for Njj = 4 
(Njj = 2) if the two polar (hydrophobic) sides are op- 
posite to each other, n c ,j = 12 for Njj — 4 (Njj = 2) if 
they are adjacent, n c j = 8 for Njj = 3 if the three polar 
sides are are all adjacent to each other, and n c j = 12 for 
Nh = 3 if two of them are opposite to each other. 

Because a water side may be in one of q different states, 
Si :t7os is 1 if it is in one of the q os ordered shell states and 
otherwise, and 5% ads is 1 if it is in one of the q ds disordered 
shell states and otherwise. Analogous considerations 
apply for the bulk states. 

We note that for completely hydrophobic solute parti- 
cles (Np = 0) Eq. CJ is not equivalent to the Hamilto- 
nian in Ref. [22] , because the sites are treated differently. 
In Eq. (J) each face of a water site i contributes to the 
free energy, and a sum is performed over all pairs of faces 
in contact with each other, as opposed to over individ- 
ual sites. The representation in terms of cube sides is 
adopted here for consistent description of lipid molecules 
in Sec. Ill Cl However, we will find in Sec. IIII Al that the 
qualitative differences between the two Hamiltonians for 
fully hydrophobic solute particles are small. 

To determine the canonical partition function, a sum is 
performed over the state configurations {<Ji}. By taking 
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into account the possible orientations of the amphiphilic 
particles through the variable k, the canonical partition 
function of the system of N sites may be expressed as 

{m} *=1 fe=l 

(2) 

where Z a = q oa e 13 E °" + qd^e ^ Ed ° for the shell (<r = s) 
and bulk (a = b) states of pure water sides. The grand 
canonical partition function of the system for variable 
solute particle number is then 

S = e " Nw Z N = e-^f^\ (3) 
N {m} 

where p represents the chemical potential associated with 
the insertion of water and N w denotes the number of 
water sites. Although the explicit terms of the model 
describe solely the states of water molecules in solution, it 
contains implicitly all multi-particle interactions between 
hydrophobic solute molecules [2^. 

The coexistence regions are characterized by measur- 
ing the UCST and the LCST for various numbers Np of 
polar sides per particle. To investigate variations in the 
effective hydrophobic interactions due to changes in po- 
larity of the solute particles, we increase systematically 
the number of polar sides per particle and determine the 
coexistence region in each case. 

The extent of aggregate formation in the system is de- 
termined from the number density of contacts between 
two hydrophobic cube sides, tih-h, between two polar 
sides, np-p, and between a polar and a hydrophobic side, 
rip-H- The number of contacts in a randomly distributed 
system with the same particle density p p is also calculated 
for comparison. In a random solution of solute particles 
with Np polar faces, whose positions and orientations are 
completely independent, these probabilities are given in 
the thermodynamic limit by 

p(H-H) = P (H) 2 , 
p(P-P) = p{Pf, 

p(P-H) = 2p(H)p(P), (4) 

where p(H) is the probability that a cube face is hy- 
drophobic, p(P) is the probability that it is polar, and 
by symmetry p(H—P) = p(P—H). In a random system, 
the probabilities of occurrence of the different faces are 
independent of the neighboring sites, and are simply 

P{H) = p p , 

z 

p(P) = (l-p p ) + ^p p . (5) 

If the contact densities tlh-h and np-p, are larger than 
their probabilities of random occurrence, and np-n, is 
correspondingly smaller, the system has formed aggre- 
gates which reduce the number of hydrophobic sides ex- 
posed to water. 



C. Lipids 

Micelles are generally formed by amphiphilic 
molecules, referred to here as lipids, which are com- 
posed of two distinct regions, the polar head and the 
hydrophobic tail. The length of the tail, which is 
typically composed of one or more hydrocarbon chains, 
is normally rather greater than the size of the polar 
head (Fig. 0. Experimental observations [11] suggest 
that lipids with tails shorter than 10 carbon atoms 
are highly soluble in aqueous solutions, while those 
whose tails exceed approximately 20 carbon atoms are 
almost completely insoluble. For those molecules with 
tail lengths in the intermediate range, which show the 
widest variety of surface-active properties, the tails are 
on average some three to five times longer than the 
dimensions of the polar head. (In lipids composed of 
two or more tails, this ratio is generally smaller.) As an 
example, the head size of the glycolipid in Figgis close 
to one third of the length of the 14-atom tail. We define 
an effective tail length I as the ratio between the tail 
length and the head size of the lipid, so that a typical 
lipid is represented by I values between two and five, 
which may also be fractional. 

For simple geometrical reasons, the total repulsion be- 
tween such a tail and the water molecules surrounding 
its sides is significantly stronger than that for the small 
tail tip of the chain. We adapt the model described in 
Sec. Ill Bl to include this aspect by assigning different en- 
ergy levels to shell water clusters in contact with the sides 
of a solute particle compared to those in contact with 
the tip. Because the number of shell water molecules 
interacting with one side of a hydrophobic tail is approx- 
imately I times that of those interacting with the tip, 
the energy associated with a site representing all of these 
water molecules is taken to be I times that for a tip site. 
A cubic solute particle in the lipid model [represented in 
Fig. HUT a)] consists of a polar head (P) and a hydropho- 
bic tail, which in turn is divided into the moderately hy- 
drophobic tip (H), situated opposite the polar head, and 
the long, strongly hydrophobic sides of the tail (S). In the 
coarse-grained model, both the tip and each long side of 
the tail are represented by a face of the cubic particle, 
but the sides interact more strongly with a neighboring 
water site than does the tip. 

The Hamiltonian of a system of N sites on a cubic 
lattice, which are occupied either by water {rii = 1) or 
by a lipid molecule (n^ =0), is then 

ff[{n;},{<^},W] = 

N 1 

i=i (ji) 

m 
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where again Xij(k) depends on the orientation state 
k of the particle at each nearest-neighbor site j of i. 
However, for the lipid model Xij(k) takes the value 
1 if the neighboring face is polar (i.e. water or the 
relevant side of a solute particle), if it represents a 
slightly hydrophobic tail tip H, and —1 if it represents 
a strongly hydrophobic tail side S. k varies again from 



1 to n c ,j, where n c j 



1 for water and n c j — 6 for 



lipid particles. Because a polar side i may be in one of 
q different states, Si !<Toa H is 1 if site i is occupied by a 
polar face in one of the q os ordered shell states of a H 
face and otherwise, and 5i^ dB H is 1 if it is occupied 
by pure water in one of the qds disordered shell states 
of a H face and otherwise. Analogous considerations 
apply for S faces and for the bulk states. We stress here 
that the Hamiltonians in Eqs. and JfJJl provide a full 
description of the microscopic states of the amphiphilic 
and lipid systems. The macroscopic, thermodynamic 
behavior is contained implicitly in these microscopic 
states [il HI El |H |p, Hi HI H H3, and the 
differences for such properties as the phase boundaries 
for varying chain length enter through the microscopic 
energy and degeneracy parameters. 



To measure the formation of aggregates in the lipid 
system, the number of contacts between different faces 
is again determined and compared with the number of 
contacts in a randomly distributed system of the same 
particle density p p . In a random solution of solute par- 
ticles, whose positions and orientations are completely 
independent, these probabilities are given in the thermo- 
dynamic limit by 



p(a—a) = p(a) 2 , 
p{a-(3) = 2p{a) P {(3) 



p(P-a), 



where a, (3 = H, P, S. Here p{H) is the probability that 
the adjacent side of the nearest neighbor is slightly hy- 
drophobic, p(P) the probability that it is polar, and p(S) 
the probability that it is strongly hydrophobic. The prob- 
abilities are given in general by 

Nu 

P{H) = Pp, 

z 

N P 

P( p ) = I 1 - Pp) H P P , 

z 

P (S) = ^ Pp: (11) 



As in Sec. Ill HI only water sites contribute to the free 
energy, and a polar face is considered having the same 
effect on a neighboring water site as water. We define 



-0E d 



Sh = --pln[q OStH e 0E °*> H + q dsH t 
S S = —H^se'^+q^se-^], 



B = -iln[g o6 e-^ + q dh e~^ 



(7) 



For a system of N sites, the canonical partition function 
may be expressed as 



z n =y, n n e^ 6 ^ x,jw(a « w+i) 

{n z } i=l fe=l 

xe -/3^[5 ff ^ Ui> (l-A? j(fc) )+ < S s ^ <ji) lA l , j(fc) (l-A 1)j(fc) )] 



(8) 



whence the grand canonical partition function for a sys- 
tem where the lipid density is not fixed is 

S = E e-W'Zs = e-^SK"*}], ( 9 ) 
N {m} 

with the effective, grand canonical Hamiltonian 

N n c ,j 



P) 



i=l k=l 



<Ji) 



although henceforth we will consider only the values 
Nh = Np = 1, Ns — 4. Aggregate formation has 
occurred in this system if tlh-h, n P~p, an d rigs are 
larger than their random expectation values, and np-n 
and rips are correspondingly smaller. 

The absolute energy levels and the degeneracies of wa- 
ter sites facing S sides are higher because this site repre- 
sents the number of water molecules contained in a shell 
site of H multiplied by the factor I. The effective energy 
of a S shell site is thus obtained from that for a H shell 
site, for both ordered and disordered states, using 



E os ,s = I E OS: h, 



E, 



ds,S 



— lEds,H, 



(12) 



where E os ,h (-E'ds.n) is the energy of an ordered (disor- 
dered) water site in the shell of the tip (H) and E os ^s 
(Ed, s ,s) that of an ordered (disordered) shell water site of 
the tail side (S). Under the assumption that neighboring 
water sites are rather independent, the total number of 
configurations for two sites may be approximated by the 
product of their numbers of configurations. Thus, the 
number of configurations of a S shell site is related to 
that of a corresponding H shell site by 



Ids.H- 



(13) 



The parameter values of the energy levels, and the de- 
generacies of bulk water and H shell sites, are chosen as 
described in Sec. Ill Al 



D. Methods 

Our interest is focused on the orientation and location 
of amphiphilic molecules in solution, which may be cap- 
tured by molecular-level simulations. We thus restrict 
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our considerations to Monte Carlo studies, because the 
processes involved depend strongly on local, spatial ef- 
fects which are neglected in mean- field calculations. As 
in Ref. [2^, we use Monte Carlo simulations to detect 
the aggregation of amphiphilic solute particles as a func- 
tion of temperature and solute concentration. We work 
with a system of 30 x 30 x 30 sites on a cubic lattice 
with random initial particle distributions and periodic 
boundary conditions, using a Metropolis algorithm for 
sampling of the configuration space. The numbers of re- 
laxation (100 000) and measurement (1 000 000) steps are 
similar to those in our previous studies j^. The closed- 
loop coexistence curves in the p p -T phase diagram are ob- 
tained from the transitions determined by increasing the 
temperature at fixed chemical potential (grand canoni- 
cal sampling), which results in a sudden density jump at 
the transition temperature, and from the corresponding 
solute particle densities. 

Analysis of the properties of a lipid-water system at 
the molecular level is possible by similar Monte Carlo 
simulations using the model of Sec. Ill CI However, the 
procedure described above must be redefined in one re- 
spect. In the Metropolis algorithm, the relative transi- 
tion probability to a configuration {n^} from a previous 
one {rii} depends on the difference in free energy of the 
two configurations according to 

r = e -/3(H£[{n' i }]--ffS[{n i }]) > ^ 

The free-energy difference must be calculated for two 
states with the same number of molecules. One step 
of the simulation procedure consists either in rotation 
of a solute particle or in a site exchange between two 
randomly chosen sites. The only contributions to the 
difference in free energy are then the energy change of 
the sites concerned, and of their nearest-neighbor sides. 
If a side changes from P to S or vice versa, the bulk 
and shell states of neighboring water sites must con- 
tain the same number of water molecules to be com- 
parable. In this case, we attribute the bulk energies 
E b,s = I Eob, Edb,s = I Edb, and the respective degen- 
eracies q b,s = 9obJ Qdb,s = Tdbi to the relevant water sites 
when calculating the probability r. The index S refers to 
the comparison of this bulk site with a S shell site. 



III. RESULTS 

A. Amphiphiles 

We begin by attempting to capture the qualitative be- 
havior of amphiphilic molecules with varying conforma- 
tions of polar and non-polar regions. Fig. [3] shows the 
Pp-T phase diagram obtained from the simple model of 
Sec. Ill Bl for different numbers of polar sides per parti- 
cle. The densities represent volume fractions, and are 
therefore dimensionless. In the coexistence regime the 
amphiphilic particles aggregate, minimizing the contact 




FIG. 3: Pp-T phase diagram for micelle formation in the 3D, 
HP model for different numbers of polar sides per solute par- 
ticle. The coexistence region is reduced, and aggregation sup- 
pressed, as the number of hydrophobic sides substituted by 
polar sides increases. At a given dimensionless temperature 
Tlcst < T < Tucst the system is homogeneous for solute 
particle concentrations below the CMC (pi), while above this 
value it separates into two phases of densities pi and p2. 



of their hydrophobic regions with water and with polar 
solute segments. Outside this region the amphiphiles are 
soluble at all densities below the LCST and above the 
UCST. For temperatures between these values, the am- 
phiphiles are soluble only at very low densities, while 
at the CMC, pi, the particles aggregate and the system 
separates into two phases: nearly pure water (of particle 
density p{) and an amphiphilic phase of density pi. With 
increasing polarity, the solubility of the solute particles 
is enhanced and the CMC increases. 

As expected, the coexistence region is reduced as the 
number of polar sides increases If the system repre- 

sents a solution of purely hydrophobic particles (Np = 
and Nh = 6), the coexistence curve is almost identical 
to that found in the HP model in Ref. J2$| (see Sec. lHBll . 
Substitution of one hydrophobic side per particle by a po- 
lar side (Np = 1 and Nh = 5) leads to a decrease of the 
UCST and a slight increase of the LCST, and to a small 
overall suppression of the temperature and density range 
of the coexistence region. In this case, which may be 
taken to represent simplified lipid molecules (discussed in 
Sec llVf) . the effect is rather moderate. However, substi- 
tution of a second hydrophobic side per particle (Np = 2 
and Nh — 4) reduces the coexistence region dramatically 
(Fig. 01 . When substituting three or more hydrophobic 
sides by polar ones, the solute particles become soluble 
at all temperatures and no aggregation phase transition 
is observed. 

Figures 0] and [5] illustrate the nature of the aggregated 
phase using "snapshots" of two-dimensional (2D) systems 
at T = 1.0, obtained in the coexistence region for so- 
lute particles with different numbers of polar sides. The 
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FIG. 4: Snapshots of 2D systems in the coexistence phase, 
obtained by Monte Carlo simulations at T = 1.0. Left: 
completely hydrophobic solute particles (Np = 0) in water; 
right: mainly hydrophobic solute particles (Np — 1) in wa- 
ter. White circles (right) represent the polar sides of the so- 
lute particles, which are shown as black squares. The parti- 
cles form compact micelles which shield the hydrophobic sides 
from water. The lattice constant a is defined by the solute 
particle size. 
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FIG. 5: Snapshots of 2D systems of partially hydrophobic 
solute particles (Np = 2) in water in the coexistence phase, 
obtained by Monte Carlo simulations at T = 1.0. 



snapshots are taken after allowing the system to relax for 
1 000 000 steps, where every 20 th step is an attempt to ex- 
change two sites and the others are attempts to rotate a 
particle. Systems containing primarily hydrophobic par- 
ticles (Np = and Np = 1) form mostly compact clus- 
ters, which minimize the number of hydrophobic surfaces 
exposed to the solvent. For lipid- like solute particles with 
Np = 1 [Fig. EI a)], the formation of perfect micelles with 
a hydrophobic core and a polar surface is prevented by 
the nature of the square lattice, which causes frustration 
on the edges of the micelles: an edge particle is forced to 
expose one of its hydrophobic sides to water. The model 
allows a lipid particle to occur in the core of micelles, 
where its polar side is in direct contact with the polar 
side of another lipid particle, because no distinction is 
made between a group of water molecules and the polar 
face of a particle. Incorporating this distinction into the 
description could be expected to generate more realistic 
micellar structures, albeit within the confines of the cubic 
geometry. 

Amphiphilic solute particles whose surface is half polar 



FIG. 6: Illustration of schematic analogs obtained using 
square particles with different arrangements of polar sides in 
2D (see text). Hydrophobic sides of the square solute sites 
are shown as thick solid lines, polar sides as dashed lines. 



(Np = Nh = 2 in 2D) show differing behavior depend- 
ing on the polarity pattern of the sides. If the polar sides 
are adjacent on the square, small micelles consisting of 
four solute particles can be formed, which is energetically 
the most favorable configuration because no hydropho- 
bic sides are exposed to water (Fig. EJ. Short, diagonal 
lines of molecules, which may be considered to represent 
condensed bilayers, can also be formed, although their 
ends are hydrophobic, and this configuration is therefore 
less favorable than are "circular" micelles of four solute 
particles. These configurations are expected from the 
construction of the sites, which is shown in Fig. El to ap- 
pear as the ground states. Solute sites with two adjacent 
polar sides may be considered to represent sections of 
circular micelles [Fig. Elk)], with the formation of small 
micelles as a consequence. In contrast, if the polar sides 
are opposite each other, the only possibility to avoid hy- 
drophobic sides being in contact with water is to form 
lines of particles, although again the hydrophobic ends 
remain exposed to water. These squares may be taken 
as schematic representations of cross-sections of a bilayer 
[Fig. E(c)]. 

The contact densities uh-h, np-p, and np_#, are 
shown in Fig. for a system of solute particles with 
one polar side (Np = 1) as a function of relative tem- 
perature. The results are normalized to the probability 
of these contacts in a randomly distributed system with 
the same particle density, p p = 0.24. At low temper- 
atures, the solute particles are clearly soluble, because 
n H -H < p(H-H), np-p < p(P-P), and n P - H > p(P-H). 
This reflects the formation of strongly hydrogen-bonded, 
partial cage- like structures of water molecules around the 
hydrophobic parts of the amphiphilic particles when en- 
tropy effects are minor. At temperatures higher than the 
lower critical temperature for density p p , increasing en- 
tropy effects favor a screening of hydrophobic faces from 
polar ones, and the solute particles aggregate to form mi- 
celles. In this regime uh-h > p(H—H), np-p > p(P—P), 
and the density of polar-hydrophobic contacts is sup- 
pressed, np^H < p{P—H). The effect on uh-h /p(H—H) is 
more pronounced than np-p /p(P—P) because the number 
of polar sides in the system is higher than the number of 
hydrophobic ones. Changes in tih-h then lead to a larger 
relative effect, and in fact the majority of P-P contacts 
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FIG. 7: Contact densities between polar and hydrophobic 
sides as a function of dimensionless temperature for a 3D sys- 
tem of solute particles with one polar side, Np = 1. The 
contact densities between sides are normalized to the values 
expected in a random system for p v = 0.24. 



are intact even in the dissolved phase due to the high 
number of water molecules. For low particle densities, 
the relative effect is larger still because fewer H-H con- 
tacts are possible. 

Finally, at high temperatures, the entropy becomes 
dominant and the contact densities approach their re- 
spective random values as complete mixing is obtained. 
We note that at high temperatures the contact densities 
do not recross the value 1 to recover the low-temperature 
phase of single-particle dissolution (Fig. 0) . Instead their 
values simply converge to unity, implying that the com- 
plete miscibility takes the form of a truly random parti- 
cle/water distribution. 



B. Lipids 

We have performed Monte Carlo simulations for lipid 
molecules in water on a cubic lattice for varying lengths 
I of the hydrophobic tail, as described in Sees. Ill CI 
and III Dl to investigate the stability of their aggregation 
as a function of hydrophobicity (I) and density. The p p - 
T phase diagram (Fig. |5J) shows clearly that the coexis- 
tence regime is enhanced significantly with increasing tail 
length /, which can be attributed to the stronger effec- 
tive repulsive interaction between the longer hydrophobic 
tail of the lipid molecule and the surrounding water. The 
ratio I can be taken as the important quantity to char- 
acterize the degree of hydrophobicity of a lipid molecule. 
Comparison with Fig. ^ indicates that for a typical polar 
head, I = 1 corresponds to a tail containing approxi- 
mately four carbon atoms. Fig. [S] shows that for lipid 
molecules with a tail containing approximately 8 carbon 
atoms (I — 2) there is already an enhancement of the 




FIG. 8: Pp-T phase diagram illustrating micelle formation in 
the extended 3D, HP model for lipid molecules of varying tail 
length. / represents the relative length of the hydrophobic tail 
compared to the size of the head (Fig. 0. The coexistence 
region is enhanced as the length of the hydrophobic tails in- 
creases, and aggregation is promoted. 



coexistence region. The aggregation of solute particles 
with a tail composed of approximately 12 carbon atoms 
(/ = 3) is reinforced very significantly. Within our sim- 
plified model, at relative temperature T = 1 lipids with 
I > 4 are basically insoluble, forming a completely sep- 
arated phase, and thus perfect micelles, for all densities, 
in rather good qualitative agreement with expectations 
based on experiments |38| . 

The aggregation may be analyzed in more detail by 
studying the density of side contacts as a function of 
temperature. The thermal evolution of the density of 
contacts between the different sides in a system of lipid 
particles of tail length / = 3 in water is shown in Fig. [5] 
The results are normalized to the corresponding prob- 
ability of these contacts in a randomly distributed sys- 
tem with the same particle density, p p = 0.25. At low 
temperatures ur-h < p(H—H), np-p < p(P—P), and 
ns-s < P(SS), while n P - H > p(P-H), n P _ s > p(P-S), 
and Tins > p(H—S), meaning that the solute particles 
are clearly dissolved as a consequence of partial cage for- 
mation around the hydrophobic tails. At temperatures 
higher than the lower critical temperature for density p p 
(Fig. [SJ , the lipid particles aggregate to minimize their 
total exposed surface, whence uh-h > p(H—H), np-p > 
p(P-P), and n S -s > p(S~S) : while n P - H < p{P-H) 
and np-s < p(P—S). Because of the much stronger ef- 
fect of an S face than of an H face on its neighboring 
water site, the normalized contact density of two S faces 
is highest in the aggregation phase and lowest below the 
LCST. A higher relative contact density is observed be- 
tween S and H faces than between two H faces for the 
same reason, while below the LCST water sites prefer to 
form the solvation shell of S faces rather than of H faces. 
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FIG. 9: Contact densities between the different faces as a 
function of dimensionless temperature for a 3D solution of 
lipid molecules with effective tail length I — 3. The contact 
densities between sides are normalized to the values expected 
in a random system for p p — 0.25. 
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FIG. 10: Representation of different surface patterns on cubic 
solute particles in the HP model, and their schematic corre- 
spondence to different micelle types(see text). Polar surfaces 
and micelle segments are marked in gray, hydrophobic sur- 
faces in white. 



At high temperatures the contact densities converge, in 
fact rather abruptly (Fig. |5J) to their random values, in- 
dicating complete mixing. 

IV. DISCUSSION 

Amphiphilic molecules in aqueous solutions can form 
different types of micelles depending on their concentra- 
tion and on the distribution of polar regions at their 
surfaces. Our initial investigation of the qualitative 
properties of micelle formation in a hydrophobic-polar 
model involved systematic substitution of the hydropho- 
bic sides of cubic solute particles by polar ones. We 
determined the p p -T phase diagram for different sur- 
face patterns and found closed-loop coexistence curves 
(Fig. [SI, in accord with experiments using hydrophilically 
modified copolymers of poly(N-isopropylacrylamide) p|. 
With increasing polarity, the coexistence region is re- 
duced as the solubility of the model amphiphiles in- 
creases. This tendency was confirmed by the same exper- 
iment, where the LCST of purely hydrophobic poly(N- 
isopropylacrylamide) P3 was observed to increase from 
37° C at atmospheric pressure to 42-44° C when approxi- 
mately 13% of the monomers were substituted by polar 
species (CP2 and CP3). In our model we observe the 
same quantitative increase of 2-3% in absolute tempera- 
ture from Tlcst = 0.545 for purely hydrophobic solute 
particles to Tlcst = 0.56 for amphiphiles with one polar 
side, which represents 17% of the particle surface. Thus 
the crude cubic model appears to yield good agreement 
with available data at this level of comparison. 

For particles with two polar sides the coexistence re- 
gion is reduced dramatically. This is not surprising, 
considering the fact that two polar sides represent one 
third of the total particle surface, and the attractive in- 



teractions with water are rather strong. In amphiphilic 
molecules the polar region is usually rather small com- 
pared with the hydrophobic surface area. The solubility 
of the molecules thus increases considerably on substitu- 
tion with polar monomers, leading to a decrease in UCST 
and an increase in LCST and CMC. In the cubic model no 
significant difference in the size of the coexistence region 
is observed for different distributions of the two sides on 
the solute cubes. Any further substitution of hydropho- 
bic faces by polar ones results in molecules which are at 
least half polar: the solubility of such particles is always 
high, and thus no aggregation is found. 

Possible schematic interpretations of the various sur- 
face patterns of cubes representing solute particles in the 
HP model are shown in Fig. The model is applica- 
ble for any surface pattern and density of solute, with 
the premise that one site may contain one amphiphilic 
molecule or a group of solute molecules. A given surface 
distribution of polar groups on a polymer may be char- 
acterized by a corresponding arrangement of polar sides 
on the surface of a cubic solute molecule. A single lipid 
molecule with a clear distinction between a polar head 
and a hydrophobic tail may be represented by a cube 
with one polar face [Np — 1, Fig. HOf al]. A small sec- 
tion of a cylindrical micelle would be represented by a 
cube with two adjacent polar sides [Np = 2, Fig. HUf c)]. 
while a particle with two opposite polar sides corresponds 
at the same level of approximation to a cross-section of 
a bilayer [Fig. HOf bl]. For the formation of "spherical" 
micelles, each site must have three adjacent polar sides 
[Np = 3, Fig. HOT d)]: if two of the three sides are situated 
opposite each other [Fig. HOf e)]. the site may again be 
considered as a section of a cylindrical micelle. Although 
by construction Np — 3 should give small spherical or 
cylindrical micelles, depending on the distribution pat- 
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tern at the surface of each site, in fact the solubility is 
too high to find aggregates. 

Micellar structures occurring in 3D simulations are dif- 
ficult to display. To confirm the formation of different mi- 
celle types depending on the surface pattern of the solute 
particles, we have considered snapshots of an analogous 
2D system (Figs. 0] and EJ). Here the square solute par- 
ticles may be interpreted in a manner similar to the 3D 
case (Fig. |HJ|, and the formation of small micelles and 
layers is found in the coexistence region. During the 
relaxation process, micelles grow from initial dimers to 
larger entities. Although the solute particles may rotate 
at a given position, they can be trapped in a configura- 
tion which disables the construction of perfect micelles 
or extended layers. Because there is no preference for 
growing a layer in one direction rather than in the other, 
short line segments are formed which are incompatible 
with others, resulting in a network of short layers. In the 
model, no distinction is made between the polar side of a 
water molecule and that of a solute particle, and a P-P 
contact contributes the same energy independent of the 
molecules to which the sides under consideration belong. 
Such contacts between the polar sides of solute molecules 
are found in the interior of a micelle, which also influences 
the formation of perfect micelles. The non-zero temper- 
ature in the coexistence region, and the observation that 
the upper critical density is much smaller than unity, 
might further imply the formation of imperfect micelles. 
In fact the extent to which the upper critical density is 
significant remains unclear, because shell water sites in 
the model may be considered as belonging to the micelle 
phase rather than to the pure water phase, which would 
explain the low density of the organic phase even for per- 
fect micelles containing no water molecules. 

We have extended our analysis to describe lipids, which 
represent a particular type of amphiphilic molecule. 
Lipids are distinguished by a special partition of the polar 
and hydrophobic segments along the molecule: a typical 
lipid molecule consists of a polar head and one or two 
hydrophobic tails. To incorporate these geometrical fea- 
tures in the model we have adapted the energy levels and 
their degeneracies according to the tail length I, where a 
typical lipid molecules exhibiting surface-active proper- 
ties would be represented by values 2 < I < 5 (38| . 

As expected from experiments || , Monte Carlo simula- 
tions of lipids of increasing length illustrate a significant 
enhancement of the coexistence region (Fig. 0. Lipids 
with a longer tail have more pronounced characteris- 
tics of hydrophobic solute particles than do lipids where 
the polar head represents a considerable fraction of the 
molecule. The solubility of long-tailed lipids is therefore 
lower than that of short-tailed ones, causing a decrease in 
LCST, an increase in UCST, and a decreased CMC at any 
given temperature as the tail becomes longer. For the pa- 
rameters used in the model we find that an increase from 
/ = 1 to I — 2 already enhances the coexistence region 
substantially, and for that I > 4 the system is completely 
separated into pure water and essentially perfect micelles 



of density p p = 1 over a wide temperature range, which 
is fully consistent with experiment |3S| . 

A quantitative measure of micelle formation in am- 
phiphilic and lipid systems is provided by the contact 
densities and n a -p, normalized by the corresponding 
probabilities of the contacts in a random system (Figs.[7| 
H5J. We have found fewer contacts between hydrophobic 
faces and more hydrophobic-polar contacts than would 
be expected for a random distribution below the LCST, 
indicating highly dissolved solute molecules. Above the 
LCST this picture is inverted, thus confirming the aggre- 
gation of amphiphiles and lipids, and the separation of 
the solution into two phases. From a knowledge of the 
system geometry and densities, the values of the contact 
ratios may also be used to confirm the extent of solu- 
tion or aggregation, and also the effective "purity" of 
the micelles which form for different polar distributions 
on the solute molecules. At temperatures in excess of the 
UCST, the solution approaches a random mixture due to 
dominant entropy effects; in Monte Carlo studies the rela- 
tive contact densities approach unity continuously, rather 
than undergoing a sharp transition. 

We close by emphasizing again the limits of our analy- 
sis. We have formulated a model for an aqueous solution 
using a minimal set of assumptions; we have used an ex- 
tremely crude representation of amphiphilic molecules, 
and have considered only a cubic system as the foun- 
dation on which our "micelles" are constrained to form. 
Nevertheless, we have obtained a realistic set of aggre- 
gation phenomena and a surprising degree of agreement 
with the available experimental results. There are, how- 
ever, several examples of phenomena which are beyond 
the reach of the model at its current level of refinement. 
A solution of lipids in water may produce a lamellar 
phase of bilaycrs at low temperature and rather high 
densities ^5|. Bilayer formation requires the possibil- 
ity of smooth curvature and high flexibility p|, and is 
precluded in the Monte Carlo simulations by geometri- 
cal constraints presented by the lattice. For the same 
geometrical reasons it is also difficult to find well-formed 
micelles in the more general case of amphiphiles in water. 

A further limitation is that the model considers ex- 
plicitly the energy states of water sites, and therefore no 
distinction is made among P-P, H-H, S-H, and S-S con- 
tacts. Because there is also no difference between a group 
of water molecules and the polar face of a solute particle, 
the model a Hows amphiphiles to occur in the core of mi- 
celles, where their polar sides are in direct contact with 
those of other amphiphiles. For lipids, an S-H contact 
which is formed during the Monte Carlo simulation is as 
favorable as an S-S contact, although it is more likely to 
be broken in a later step. In addition, the orientation of 
a lipid molecule is irrelevant for the formation of an S-S 
contact, which prevents efficient alignment of the heads 
and may further hinder the formation of extended bilay- 
ers. The incorporation of these distinctions in a more 
sophisticated description of the lipid solution may be ex- 
pected to reproduce further detailed properties of real 
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systems. 

In summary, we have extended the MLG framework 
to include the solvation of amphiphilic solutes in water. 
Within a cubic HP model we have found the aggregation 
of solute particles, and the formation of various types of 
micelle as a function of the distribution of hydrophobic 
regions. By successive substitution of hydrophobic sides 
by polar ones, we have studied the aggregation behavior 
and the influence of the degree of hydrophobicity on 
the upper and lower critical solution temperatures. 
We have refined this model to describe lipid molecules 
of varying length, by adapting the interaction of the 
hydrophobic tail to include a corresponding number of 
neighboring water molecules, and have demonstrated 
the enhanced stability of aggregates with increasing 
tail length (increasing hydrophobicity). We have shown 



that primary features of micelle formation, which are 
often attributed solely to the amphiphilic nature of the 
solute particles under consideration, are reproduced 
by our extension of the solvent-based MLG model to 
describe alterations of water structure in the vicinity of 
the different surface regions of dissolved amphiphiles. 
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